Features due to spin-orbit coupling in the optical conductivity of single-layer graphene 
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We have calculated the optical conductivity of a disorder-free single graphene sheet in the presence 
of spin-orbit coupling, using the Kubo formalism. Both intrinsic and structural-inversion-asymmetry 
induced types of spin splitting are considered within a low-energy continuum theory. Analytical 
results are obtained that allow us to identify distinct features arising from spin-orbit couplings. We 
point out how optical-conductivity measurements could offer a way to determine the strengths of 
spin splitting due to various origins in graphene. 
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I. INTRODUCTION 

Graphene is a single sheet of carbon atoms form- 
ing a two-dimensional honeycomb lattice. This mate- 
rial has only recently become available for experimen- 
tal study, and its exotic physical properties have spurred 
a lot of interest.—!^ Known theoretically since the late 
40s,— graphene is a promising candidate for applica- 
tions due to its excellent mechanical properties;^ scal- 
ability down to nanometer sizes, ^ and exceptional elec- 
tronic propertiesi^ The conical shape of conduction and 
valence bands near the K and K' points in the Bril- 
louin zone renders graphene an interesting type of quasi- 
relativistic condensed-matter systeni^'^ where mass-less 
Dirac-fermion-like quasiparticles are present at low en- 
ergy. In contrast to the truly relativistic case, the spin 
degree of freedom in their Dirac equation corresponds 
to a pseudo-spin that distinguishes degenerate states on 
two sublattices formed by two nonequivalent atom sites 
present in the unit cell. 

The pseudo-spin degeneracy can be broken by spin- 
orbit interaction (SOI), which mixes pseudospin and 
real spin. There has been huge interest in SOI 
in graphene, resulting in a large body of theoreti- 
^,^19^^102.13^05.1607^8^^^,22^^ and experimen- 

tal^i^'^i^'^ work. There are two main causes for 
the SOI in graphene. Firstly, external electric fields 
(e.g., due to the presence of a substrate, a backgate, 
or adatoms) and local curvature fields (ripples) induce 
a SOIi2iO^Oii2£.2i whose coupling strength we denote 
by Afj. We refer to this contribution as the Rashba 
SOI in the following. In addition, there is an intrinsic 
SOia02030i05|i6i24 ^^^-^ strength A/, which is caused by 

the atomic Coulomb potentials. 

Existence of the intrinsic and Rashba SOIs can be in- 
ferred from group-theoretical arguments i^O^i^S However, 
the actual values of their respective strengths and 
A/ are the subject of recent debate. Initial estimates^^ 
have been refined using tight-binding models^^'"'^'* and 
density-functional calculations O^OSii^ First experimen- 



tal observations of spin-orbit-related effects in graphene's 
band structure based on ARPES data^^i^ have later 
been questioned. ^^'^^ Detailed knowledge about typical 
magnitudes and ways to influence A/j and A/ is crucial, 
e.g., for understanding spin-dependent transport^ and 
spin-based quantum devices^ in graphene. The desire 
to identify possible alternative means of observing, and 
measuring, spin-orbit coupling strengths in graphene has 
provided the motivation for our work reported here. 



We present a theoretical analysis of graphene's 
optical conductivity cr(w), extending previous stud- 
jgg32^^^^^^^^ ^jjg situation with finite SOI. 

SOI effects on the DC conductivity were investigated in 
a recent theoretical study for a bipolar graphene pn junc- 
tion^ and the effect of intrinsic SOI on the polarisation- 
dependent optical absorption of graphene was considered 
in Ref. J^. Our study presents the analogous scenario for 
the richer case of the optical conductivity when both in- 
trinsic and extrinsic types of SOI are present. Since A^ 
can be tuned by external fields, we will analyze various 
situations distinguished by the relative strengths of An 
and Aj. 



Our findings suggest that optical-conductivity mea- 
surements can be useful to identify and separate different 
SOI sources. We work on the simplest theory level (linear 
response theory, no interactions, no disorder) and disre- 
gard boundary effects for the moment. The structure of 
the remainder of this article is as follows. In Sec. |TT1 we 
summarize basics of our calculation of the optical conduc- 
tivity based on the Kubo formalism; except for some de- 
tails that have been relegated to an Appendix. In Sec. lIIIl 
we show results for different relative magnitudes of SOI 
strengths at finite temperature T and chemical potential 
/i. Finally, in Sec. IIVI we summarize our results and 
discuss their applicability to actual experiments. 
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II. OPTICAL CONDUCTIVITY 



where a,b = x.y and the kernel reads 



We start from a low-energy continuum description of 
graphene^ H{k) = Hoik) + Hr + Hi. Without the SOI 
terms Hr and if/, the single-particle Hamiltonian in 
plane-wave representation reads 



(1) 



with Fermi velocity v ~ 10® m/s. The Pauli matrices a^^y 
act in pseudo-spin space, where the two eigenspinors of 
az correspond to quasiparticle states localized on sites 
of the A and B sublattice. Analogous Pauli matrices 
Tx,y,z act in the two-valley space spanned by states near 
the two K points. The part of the effective Hamiltonian 
describing Rashba SOI is given by 



J^R = yCTxTzSy — (TySx), 



(2) 



where Ar includes both the external electric-field and 
curvature effects in a coarse-grained approximation, with 
the latter assumed to be homogeneous. The Pauli matri- 
ces Sx,y,z act in the real spin space. For the intrinsic SOI 
induced by atomic potentials, we have 



Hi = AiOzTzSz 



(3) 



The full Hamiltonian is then an 8 x 8 matrix in the com- 
bined sublattice, spin, and valley space. 

The full Hamiltonian matrix turns out to be block- 
diagonal in the valley degree of freedom, and each block 
can be transformed into the other via a unitary trans- 
formation. The bulk spectrum - ignoring subtleties re- 
lated to the topological insulator phase encountered for 
2A/ > Afl/ - for now - can then be obtained from a 4 x 4 
Hamiltonian matrix in the basis (A 1,5 j) at 

one K point. The valley degree of freedom then merely 
manifests itself as a degeneracy factor — 2. The energy 
spectrum is obtained as 

£k,..' = \ [i^'Ar + iy^4{hvy I k |2 +{Ar - 2v'Ai) 

(4) 

where the combined indices v^v' = ±1 label the four 
bands. The corresponding eigenstates 



= |k) 



(5) 



are composed of a plane wave state |k) and a k-dependent 
4-spinor \l'v')■^^■ 

We compute the optical conductivity using the stan- 
dard Kubo formula4^ 



le 

Kab = 1-Tr 

h 



(7) 

Here e denotes the electron charge, rf, is a Cartesian com- 
ponent of the position operator, po the equilibrium den- 
sity matrix, and the current operators are given by 



i e dH{k) 
Ja^ ^ [H{k),ra] - 



h dka 



(8) 



Following Ref. |32, we use the single-particle eigenstates 
\n) and eigenenergies e„. The conductivity then reads 



{n\[H,ra]\n') {n'\[H,n]\n) 



i ^ (e„/ + fiu) - zO+) 



X [/(£„) _/(e„,)] 



(9) 



where /(e) is the Fermi-Dirac distribution containing the 
chemical potential [i and the inverse temperature — 
l/iksT). 

In the absence of a magnetic field, the off-diagonal en- 
tries vanish, axy — 0, while symmetry arguments show 
that axx = cTyy = cr(uj). At finite uj in the clean system, 
only the inter-band contribution to the conductivity is 
relevant. Its real part is given by 



Re a{uj) — ire 



X 



(27r)2 

/(£k,Ki/) ~ /(ek.Kv) 



X [(5(ek,Ki/ - ek.K'i/' + fii^) + ^(ekKV - £kKi/ + ftw)] 



(10) 



where 



<.V.'(k)=k(A^i^ba|«V')i 



are the current operator matrix elements in the eigenba- 
sis, and Y! = ll(^vmn'u')- We also used 

<,,,,,,(k) = K,,,^,,(k)]^ 

since the current operator is Hermitian. In what follows, 
we restrict ourselves to the real part of (j{lo) and omit 
the "Re" sign. 

The result obtained for a; > can be expressed very 
generally as 



-2^^i^„(cc;,A,j,A/,/3,/i), 



(11) 



n=l 



<Jab 



dt e^("-'°^)*if„fc(t), 



where uq = g^e^ / (2^11) and, with the Heaviside function 
(6) O, the quantities Fn are given by 
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Fi 
F2 

Fz 

F4 

F5 
Fe 



Fi Q{hLu~ |Afl-2A/|) 

F2 [e(Afl - 2Ai)e{hLu - AR)e{2Ai + Ar - hiu) + e(2A/ - AB.)eihuj - Afl)e(2Afl - huj) 
F3 \e{AR - 2A/) + e(2A/ - AR)Q{huj - 2A/ + A^, 
F4 Q{Ar - 2A/)e(fiw - 2Ar) + e(2A/ - AR)e{hLo - 2Ai - Ar) 
F5 Q{Ar - 2Ai)e{hLu ~Ar + 2A/)e(Afl ~hLu) + e(2A/ - Ar)Q{Ar - hw) 
Fg Q{nw-2Ai- Ar). 



(12) 



r 



h 












1-3, 5-7, 10, 11 


6-8 


1-3, 5-7, 9-11 


1, 2 


10-12 


1, 2, 5 



TABLE I: List of Fn functions and the regions in which they 
contribute, as illustrated in Fig. [T] 
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FIG. 1: Regions in the An/huj — 2A///ia;-plane where the 
different contribute, cf. Table U There is no contribution 
to (7{uj) from region 4. 



The rather lengthy analytical expressions for the quan- 
tities F„{u}, Ar, A/, P, /i) can be found in the Appendix. 
In Fig. [1] we show the regions in the Ar/Huj — 2Aj/tiuj- 
plane where the different F„ contribute. 



III. RESULTS 

We now discuss the main physical observations aris- 
ing from Eqs. (jlip and (dH). First, the behavior of the 
conductivity is qualitatively different in the two regimes 
Ar > 2 A/ and Ar < 2A/. It is well-known that the lat- 
ter regime corresponds to a topological insulator phase 
while the former yields a conventional band insulator, 
with a quantum phase transition in between. For the 



topological insulator phase 



12.44.45.46 



spin-polarized gap- 



less edge states forming a helical liquid will dominate the 
optical conductivity when both fcsT and fku are smaller 
than the gap energy. In that regime, the conductivity 
is expected— to exhibit power-law behavior analogous 
to that found for ordinary one-dimensional electron sys- 
tems 4^1^ In what follows, we consider the frequency and 
temperature range such that the optical conductivity is 
still mostly determined by the bulk states. 

Sharp features are exhibited by the conductivity as a 
function of frequency uj, which depend on the relative 
strength of the two SOI terms and should therefore allow 
for a clear identification of these couplings. We start by 
discussing a few special cases. For Ar = but finite 
A/, the gapped spectrum consisting of two doubly (spin- 
)degenerate dispersion branches leads to a vanishing con- 
ductivity for huj < A/, and all other features expected in 
the presence of a generic mass gapi^i^ In contrast, for 
A/ = but finite A^, the band structure mimics that 
of bilayer graphene, only with a gap smaller by up to 4 
orders of magnitudei^i^ The optical conductivity for this 
case has the same functional form as the conductivity for 
bilayer graphene, except that the McClure^^ inter- 
layer hopping constant is replaced by A^;. In particular, 
it exhibits a (5-peak at huj = A^j and a kink at huj — 2Ar. 
With hvk — e, the analytical expression is 



a 



-S{hw - Ai 



de 



eA, 



4e2 



A, 



4e2 



A^ 
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Ar 



'huj + 2AR haj-2Aj 



huj 



hw 



■Q{huj-2AR) 



l huJ + Ar 

\ 2 



4e2 + Al 

huj - Ar 



(13) 
Q{niu-AR), 
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FIG. 2: Optical conductivity at T = 1 K for graphene with 
A_R — 100 /icV and A/ — 0.2A_r, thus realizing the case 
2A/ < Aij. We set fi = A/ to maintain charge neutrality. 
Inset: Low-energy part of the bandstructure. Kinks in the fre- 
quency dependence of cr arise when new transitions between 
different bands become possible at certain critical values of 
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FIG. 3: Optical conductivity at T = 1 K for graphene ob- 
tained for the special case 2A/ = A_r with A_r — 100 /leV, 
setting /i — Ai = Ah/2 to ensure charge neutrality. Inset: 
The bandstructure shows that three bands cross at = and, 
hence, some of the kinks present in Fig. [2] disappear. 



where we define the function 



cosh(/i/3) + cosh(e/3) ' 

In the hmit A^j = A/ = 0, the optical conductivity of 
clean graphene with its spin-degenerate linear dispersion 
is recovered. '^^i^ The asymptotic behavior for large fre- 
quencies turns out to be independent of the SOI cou- 
plings, with a always approaching the well-known uni- 
versal value e^/(4/i). 

The optical conductivity for various situations where 
both Aij and A/ are finite is shown next in a series of 
figures. In particular. Fig. [5] shows the case where An > 
2Aj. In Fig. [3l we are at the special point A/j = 2Aj. 
Furthermore, Fig. d] illustrates the regime where A^ < 
2A/. To be specific, all these figures are for T = 1 K. 
Finally, Fig. [5] displays the effects of thermal smearing. 

For 2Aj < Aj^, we observe a splitting and widening 
of the (5-peak at A/j, while the kink at 2Afl stays at the 
same position. In addition, we observe kinks at Afl±2A/, 
see Fig. [2j At the quantum phase transition point 2 A/ = 
Afl, the dispersion exhibits a crossing of two massless 
branches with a massive branch, see inset of Fig. [3] As 
a consequence, certain sharp features exhibited by the 
optical conductivity in other cases disappear . For 2A/ > 
An, see Fig.m the conductivity shows kinks at hut — An, 
At huj = 2An, and at 2A/ ± A^;. 

We have chosen to show a very wide range of SOI 
parameters A^ and A/ in these figures. Previous es- 
timates for these parameter o-'^^'^'^i^^'^^i"'^^'^"^ range from 
0.5 ^leV to 100 ^ieV for A/, and 0.04 //eV to 23 ^ieV 
for An- The Rashba coupling is expected to be linear 
in the electric backgate field, with proportionality con- 
stant 10 /zeV nm/V (Ref. I24h . allowing for an experi- 
mental lever to sweep through a wide parameter range. 



On the experimental side, the picture is currently mixed. 
One recent experimental study^^ finds A^^ = 370 /leV 
(210 /icV) for electrons (holes) in carbon nanotubes. A 
much larger value A^ = 13 meV has been reported for 
graphene sheets fabricated on a nickel surface. 

For low temperatures (e.g., at T = 1 K in the above 
figures), the SOI couplings can be distinguished by the 
different peak structures appearing in the optical con- 
ductivity. Increasing the temperature leads to thermal 
smearing of those features, as illustrated in Fig. [5l How- 
ever, the characteristic SOI-induced peak and kink fea- 
tures should still be visible in the optical conductivity up 
to r « 10 K, albeit with a smaller magnitude. 



0.8 




huj/AE 

FIG. 4: Same as Fig. [2] except that A/ = Q.^Ar with A^ = 
100 /ieV, thus realizing the case 2 A/ > A_r. Charge neutrality 
is maintained by setting /i = Ah/2. 
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FIG. 5: Same as Fig.[5J focusing on the region 0.4 < hcj/ Au < 
1.6. The solid curve is for T — 1 K, the dashed curve for 
T = 10 K, and the dot-dashed one for T = 100 K. The dis- 
tinct kinks are thermally smeared and suppressed at elevated 
temperatures, but remain visible up to T « 10 K. 



IV. CONCLUSIONS 

We have calculated the optical conductivity for a 
graphene monolayer including the two most relevant 
spin-orbit couplings, namely the intrinsic atomic contri- 
bution A/ and the curvature- and electric-ficld-induced 
Rashba term Aj^. Our result for the optical conductivity, 
which we presented for finite temperature and chemical 
potential, shows kinks and/or peaks at frequencies cor- 
responding to Aji, 2 All, and |Ajj±2A/|. Measuring the 
optical conductivity in a frequency range covering these 
energy scales can be expected to yield detailed insights 
into the nature of spin-orbit interactions in graphene. 

We did not analyze disorder effects but expect all sharp 
features to broaden since the (5-functions in Eq. ^TU\i effec- 
tively become Lorentzian peaks. We also did not consider 
the effect of electron-electron interactions. While renor- 
malization group studies indicate that weak unscreened 
interactions are marginally irrelevant interactions may 
still play an important role. For instance, Ref. ,50, consid- 



ers interaction effects on the optical properties of doped 
graphene without spin-orbit coupling. Interactions cause 
inter-band (optical) and intra-band (Drude) transitions 
and thus a finite DC conductivity. We expect that the 
peak and kink structures arising from the spin-orbit cou- 
plings survive, however, because the relevant contribu- 
tions are additive. 

Recent experimental studies suggest that an optical 
measurement of the conductivity in the energy range rel- 
evant for SOI should be possible. Fei et al^ have mea- 
sured the optical conductivity from hw = 1.54 eV up 
to 4.13 eV. Sfightly lower energies (0.2 eV to 1.2 eV) 
were reached in Ref. [H. We suggest to perform low- 
temperature experiments at microwave frequencies, with 
energies ranging from several /leV to a few meV. 
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APPENDIX A: DEFINITION OF AUXILIARY 
FUNCTIONS 

Here we provide the six functions Fn{uj, Aji, A/,/3,/i) 
(with n — 1, . . . , 6) entering Eq. We use the follow- 

ing abbreviations: 

= i (Afl - v/(Afl - 2Az)2 + Ay^) , 

£2(2/) = I [Ar + ^{Ar - 2Az)2 + 4y2) , 

esiy) = \ (-Ar - v/(Afl + 2A,)2+4j/2) , 

= \ [-^R + v/(Afl + 2A,)2+4j/2) . 

Furthermore, we define the quantities (setting here fi, = 1 
for simplicity) 
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2/1 = 

y2 = 

2/3 = 

Vi = 

y5 = 

ye = 



i^-4A| + 4A,AH-A|, + a;^ 

^^jySA^Afl - 2A| - 4A|a; + SA^o; - 4Afla;2 + ^3 

V4A^ - SAflW + 4a;2 
V-SA^^Aflo; + 2A^a> - 4Aja;2 _^ 5^^0,2 ^ 4Afla;3 +IJ4 
2v/A^ + 2Ai^t^ + t^2 

Vw^SA^Afl - 2A^ - 4A^a) + SA^w - 4ARa;2 _,_ 



V4A2j-8A7ja) + 4a)2 
Vw^8A2 Afl _ - 4A2u; + SA^^w - 4Aflw2 ^ ^3 



V^4A2, - SAflW + 4a;2 



-4A2-4A/Afl- A|j + a;2. 



Finally, we define A± = A/j ± 2A/. With these conventions, the functions F„(a;, Ajj, A/, /3, yu) can be expressed 
follows: 



as 



A = [/(ei(j/0)-/(e2(yi))] 



= [/(ei(y2))-/(e3(2/2))] 



16 (42/2 + A^) 



#2 



3/2 



4y2 + A2 



41/2 



V4y2 + A2_ ^4y2 + A2 



2/i 


(-2A;^. + ^4y2 + A2^ - y^4y2 + A^] 




(4y2 + A_ 




-^4y2 + A2_)) (4y2 + A+ (A+ + ^ 


/42/i + A^)) 



i^3 = 



[/(ei(t/3))-/(e4(2/3))]x 

y32^4y32 + A2_ ^4^32 + A2_ (-2A« + ^4j/2 + A2_ + ^^yl + A^) 

4 (^4^/32 + A2_ + ^4^2 + A^) (4^32 + A_ (a_ - ^4y2 + A2_)) (4^2 + A+ (a+ - ^4y2 + A^)) ' 
[/(e2(2/4))-/(e3(y4))]x 

n\ sl^li\ + Al /i;;[T^ (2A,. + ^/4,2 + + ^/4/;2 + A^) 

4 (^4j/| + A2_ + ^4y2 + A2+) (4^! + A_ (a_ + ^4y2 + A2_)) (4^2 + A+ (a+ + ^4y2 + A^)) ' 

42/5 4^5 



= [/(e2(2/5)) - /(e4(2/5))] 



/4j/2 + A2_ ^4y2 + A^ 
yl (2A« + ^4y2 + A2. - ^4^/2 + A^) 



(4yf + A_ (a_ + ^4j/2 + A2_)) (4^2 + A+ (a+ - ^i^f^A^) ' 

2/6A^ 



[/(e3(2/6))-/(e4(2/6))] 



16 (4%2 + A^) 



3/2 



42/1 + A^ 



2/6 
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